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Abstract 

The lowest order constrained variational (LOCV) technique has been used to investigate some 
of the thermodynamic properties of spin polarized hot asymmetric nuclear matter, such as the free 
energy, symmetry energy, susceptibility and equation of state. We have shown that the symmetry 
energy of the nuclear matter is substantially sensitive to the value of spin polarization. Our 
calculations show that the equation of state of the polarized hot asymmetric nuclear matter is stiffer 
for the higher values of the polarization as well as the isospin asymmetry parameter. Our results 
for the free energy and susceptibility show that the spontaneous ferromagnetic phase transition 
cannot occur for hot asymmetric matter. 
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I. INTRODUCTION 



The possible occurrence of the spontaneous phase transition to a ferromagnetic state in 
nucleonic matter is very important for studies relevant to astrophysical problems, such as 
the physical origin of the magnetic field of the pulsars jl-4|. The study of the thermody- 
namic and magnetic properties of spin polarized hot asymmetric nuclear matter, such as 
free energy, symmetry energy, magnetic susceptibility plays a crucial role in understanding 
the ferromagnetic phase transition, the equation of state and the structure of systems as 
diverse as the neutron rich nuclei and protoneutron stars. 

A protoneutron star (newborn neutron star) is born within a short time just after the 
supernovae collapse. In this stage, the interior temperature of the neutron star matter is in 
the order of 20-50 MeV ^. Determination of the dependence of the symmetry energy on 
the spin polarization and the behavior of the magnetic susceptibility versus the density are 
of special interests in description of the occurrence of the ferromagnetic phase transition in 
asymmetric nuclear matter. In addition these quantities are useful to estimate the mean 
free path of the neutrino in the dense nucleonic matter which is a relevant information for 
understanding the mechanism of the supernova explosion and the cooling process of the 
neutron stars Q|. 

There may exist several possibilities for the generation of the magnetic field in a neu- 
tron star. Among them are the conservation of the magnetic flux of the original star, a 
kind of dynamo mechanism, phase transition to a ferromagnetic state, or a combination of 
above mechanisms. The possibility of the existence of a phase transition to a ferromagnetic 
state in the neutron matter and nuclear matter has been studied by several authors [7-32], 
without any general agreement. In most calculations, the neutron star matter is approxi- 
mated by pure neutron matter, at zero temperature. Some calculations such as those based 



on the hard sphere gas model [7|, 1 8| , Skyrmelike interactions [21|, variational calculation 
using the Reid soft-core potential [iJ] and relativistic Dirac-Hartree-Fock approximation 
with an effective nucleon-meson Lagrangian 19|, show that the neutron matter becomes 
ferromagnetic at some densities. There are other calculations such as Monte Carlo 23| and 



Brueckner-Hartree-Fock calculations 



24 



271 ] using modern two-body and three-body realis- 



tic interactions, which show no indication of ferromagnetic transition at any density for the 
neutron matter and asymmetric nuclear matter. Properties of the polarized neutron matter 
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32|. Bombaci et al. 32| have 



at finite temperature have been studied by several authors [30 
studied these properties within the framework of the Brueckner-Hartree-Fock formahsm us- 
ing the AVis nucleon-nucleon interaction. Their resuhs show no indication of a ferromagnetic 



transition at any density and temperature. Lopez- Val et ah |3l| have used the Dl and DIP 
parameterization of the Gogny interaction and the resuhs of their calculation show two dif- 
ferent behaviors. Whereas the DIP force exhibits a ferromagnetic transition at a density of 
p ~ whose onset increases with temperature, no sign of such a transition is found 

for Dl at any densities and temperatures. Rios et al. 30| have used Skyrme-like interactions 
and their results indicate the occurrence of a ferromagnetic phase of the neutron matter. 
The influence of the finite temperature on the antiferromagnetic (AFM) spin ordering in the 
symmetric nuclear matter with the effective Gogny interaction, within the framework of a 
Fermi liquid formalism, has been studied by Isayev {33! . 



34j . Here in our article, we use the 



lowest order constrained variational (LOCV) formalism to investigate the possibility of the 
transition to a ferromagnetic phase for the polarized hot asymmetric nuclear matter. 

he LOCV method has been developed to study the bulk properties of the quantal fluids 



35 



37|. This technique has been used for studying the ground state properties of the finite 



nuclei and treatment of isobars 



38l440l|. Modarres has extended the LOCV method to the 



finite temperature calculations and has applied it to the neutron matter, nuclear matter 



and asymmetric nuc 



of these systems 41 



ear matter in order to calculate the different thermodynamic properties 
-|44|. A few years ago, we calculated the properties of nuclear matter 



at zero and finite temperatures using the LOCV method with the new nucleon-nucleon 



potentials 
matter 



45l-l47j . Recently we have computed the properties of the spin polarized neutron 



48| , the spin polarized symmetric |49|] and asymmetric nuclear matters and neutron 



star matter j50| at zero temperature. In these works the microscopic calculations employing 
the LOCV method with the realistic nucleon-nucleon potentials have been used. We have 
concluded that the spontaneous phase transition to the ferromagnetic state does not occur. 
We have also calculated the thermodynamic properties of the spin polarized neutron matter 



5l| and symmetric nuclear matter 52] such as the free energy, magnetic susceptibility, 



entropy and pressure using the LOCV method at finite temperature. Our calculations do 
not show any transition to a ferromagnetic phase for hot neutron matter and hot symmetric 
nuclear matter. 

In the present work, we want to calculate the properties of spin polarized asymmetric 
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nuclear matter with the LOCV technique at finite temperature employing the AVis potential 
53|. 



II. LOCV CALCULATION OF THE SPIN POLARIZED HOT ASYMMETRIC 
NUCLEAR MATTER 

Spin polarized asymmetric nuclear matter is an infinite system composed of spin-up and 
spin-down neutrons with densities p^^^ and p^~\ respectively, and spin- up and spin-down 
protons with densities and p^~\ respectively. The total densities for neutrons (p„) and 
protons (pp) are given by: 

Pp = p'^^+Pi~\ Pn=pl^^+pl-\ (1) 

and the total density of the system is 

p = Pp + Pn- (2) 

Labels (+) and (-) are used for spin-up and spin-down nucleons, respectively. One can use 
the following parameter to identify a given spin polarized state of the asymmetric nuclear 
matter, 

^ PtlfiL^ S^ = p1^ (3) 

Pp Pn 

5p and 5„ are proton and neutron spin asymmetry parameters, respectively. These parame- 
ters can have values in the range of 0.0 (unpolarized) to 1.0 (fully polarized). The asymmetry 
parameter which describes the isospin asymmetry of the system is defined as, 

/3 = Pn^. (4) 

Pure neutron matter is totally asymmetric nuclear matter with /3 = 1 and the symmetric 
nuclear matter has /3 = 0. 

To obtain the macroscopic properties of this system, we should calculate the total free 
energy per nucleon, F = E — TS, where T is the temperature, and E and S are the total 
energy and entropy per nucleon, respectively. In the case of spin polarized asymmetric 
nuclear matter, the free energy per particle can be calculated by a parabolic approximation 
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resulted from the charge independence and time-reversal invariance of the nucleon-nucleon 



interaction as follows 



25|, 



F(p, r, /3, 5„, 6p) = F„{p, r, /3 = 0, Sn = 0, 5p = 0) + Flip, T){^^6n + ^-^^f 



+F,{p, r)/3^ + F,{p, r)(^5„ - (5) 

where the coefficients Fi{p, T), F2{p, T) and F^{p, T) have been determined in the following 
way, 

Fi(p, T) = F(p, r, /3 = 0, 5„ = 1, 5p = 1) - F,„(p, r, /3 = 0, 5„ = 0, 5^ = 0), 
F2(p, T) = F(p, r, /3 = 1, 5„ = 0, 5p = 0) - F,„(p, T, /3 = 0, 5„ = 0, 5^ = 0), 
F3(p, T) = F(p, r, /3 = 1, 5„ = 1, 5p = 1) - F„(p, r, /3 = 0, 5„ = 0, 5^ = 0) 

-Fi(p,r)-F2(p,r). (6) 



We calculate the total energy per nucleon {E) using LOCV method as follows |5l|, l52 |. 
We adopt a trial many-body wave function of the form 

= ^0, (7) 

where is the uncorrelated ground state wave function of A independent nucleons (simply 
the Slater determinant of the plane waves) and J-" = ^-"(1 ■ ■ ■ A) is an appropriate A-body 
correlation operator which can be replaced by a Jastrow form i.e., 

J^ = Sl[f{zj), (8) 

i>j 

in which iS is a symmetrizing operator. Now, we consider the cluster expansion of the energy 



functional up to the two-body term 



a 



For the polarized hot asymmetric nuclear matter, the one-body term Ei is 

^1 = E E Efj, (10) 

j=p,n i=+,- 

where Eij is the one-body energy of nucleon j with spin projection i, 

E^ = T.^4\k,r,f^^ (11) 

k 



nj'\k,T, pf) is the Fermi-Dirac distribution function, 



.(0 



g/3[.W{fc,r,p«)-/.«(r,pW)] 



-1 



(12) 



In the above equation P = , pf* is the chemical potential of nucleon j with spin 



projection i which is determined at any values of the temperature (T), number density (p 
and polarization {5j) by applying the following constraint, 



$:nW(fc,r,pf) = 4 

k 



(0 



(13) 



and e^*^ is the single particle energy of nucleon j with spin projection i. In our formalism, 
the single particle energy of nucleon j with momentum k and spin projection i is written 



approximately in terms of the effective mass as follows 



31 



3^ 



2m*f\p,T) 



+ u'>ir,p''). 



(14) 



p)') is the momentum independent single particle potential. In fact, we use a 
quadratic approximation for single particle potential incorporated in the single particle en- 
ergy as a momentum independent effective mass and introduce the effective masses, m*^*\ as 



variational parameters 51 



54| . We minimize the free energy with respect to the variations in 



the effective masses and then obtain the chemical potentials and the effective masses of the 
spin-up and spin-down nucleons at the minimum point of the free energy. This minimization 
is done numerically. 

The two-body energy E2 is 



E2 = ^Y.{^JHU)\^J-J^), 



(15) 



where 



2m 



[/(12),[VL,/(12)]] + /(12)\/(12)/(12). 



(16) 



In above equation, /(12) and V{12) are the two-body correlation and potential. In our 
calculations, we use the AVis two-body potential which has the following form |53l |. 

18 



Vil2) = J2v^^\n2)0 
p=l 



(p) 

12 ) 



(17) 
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where 



^(p=l-18) 



'12 



1, CTi ■ 0-2, T"! ■ T2, ((Ti • (72) (ti ■ T2) , Su, S^in " 7'2) 

L ■ S, L ■ S(ri ■ r2), L^ L2(ai ■ ^2), L2(ri ■ r^), 
L\a^■cx2)ir^■r2), (L ■ S)^ (L ■ S)2(ri ■ r2), 
T12, (o"! • (72)Ti2, 512X12, (rzi + rz2). 



(18) 



In above equation, 5*12 = [3(cri • f){a2 ■ f) — oi • a2\ is the tensor operator and T12 = 
[3(ti ■ r)(r2 ■ r) — ri ■ r2] is the isotensor operator [53|. In the LOCV formalism, the two-body 
correlation /(12) is considered to have the following form 37l |. 



/(12) 



E/^'nri2)Pi?, 
fe=i 



(19) 



where 



p(fc=l-3) 
-'12 



1 1 



1(2) 



4 4^12 ; , 



1 1 



(2) 



2 + e"- 



\0 



(5)' 
12 / , 



1 



^0(5) 



12 



(20) 



'1 + lo(^) 

.4 12 6 

The operators Oy^ and O12'' are given in Eq. f|T8|) . 

In our formalism, we minimize the two-body energy E2 with respect to the variations in 
the correlation functions fo!"^^ subject to the normalization constraint 



37 



46| 



A 



(21) 



In the case of polarized nuclear matter, the Pauli function hs^{r) is as follows 



1 _ 1 /f!M 



-1/2 



±1 





where 



7«(r) 



2z/ 
(2^ 



(22) 



(23) 



Here is the degeneracy of the system. From the minimization of the two-body cluster 
energy, we get a set of coupled and uncoupled differential equations the same as presented 
in Ref. 52|. We can obtain the correlation functions by solving these differential equations 



and then calculate the two-body energy. Finally, we can compute the energy and the free 
energy of the system. 



III. RESULTS AND DISCUSSION 



Fig. [T] shows the free energy per nucleon of the spin polarized hot asymmetric nuclear 
matter versus the total number density (p) for different values of neutron polarization 
proton polarization {6p) and isospin asymmetry parameter at T = 10 and 20 MeV. It 
can be seen from this figure that at each temperature for a given value of isospin asymmetry 
parameter, the free energy of polarized hot asymmetric nuclear matter increases by increasing 
the polarization. For all temperatures and isospin asymmetry parameters, we do not see any 
crossing between the free energy curves of different polarizations. The difference between 
the free energies of different polarizations increases by increasing the density. Therefore, we 
can conclude that there is no spontaneous transition to the ferromagnetic phase for the hot 
asymmetric nuclear matter. Fig. [T]also shows that for a given polarization and temperature, 
the free energy increases by increasing the isospin asymmetry parameter. From Fig. [1], we 
see that the free energy of spin polarized asymmetric nuclear matter decreases by increasing 
the temperature. We have found that only for the lower values of temperature (T), neutron 
polarization proton polarization (Sp) and isospin asymmetry parameter the free 
energy curve shows a minimum point. However, at higher values of these quantities, this 
minimum point disappears showing no bound state for spin polarized hot asymmetric nuclear 
matter. 

For the unpolarized case of nuclear matter with /3 = 0.0 (symmetric nuclear matter), we 
have compared the free energies at different temperatures in Fig. O It is seen that the free 
energy of unpolarized symmetric nuclear matter decreases by increasing the temperature, 
especially at low densities. From Fig. |21 we can see that at zero temperature, the saturation 
density (density of minimum point of energy) for the unpolarized symmetric nuclear matter 
is about p = 0.31 fm~^ which is greater than the empirical value (0.16 fm~^). Our results 
also show that at T = MeV , the energy of symmetric nuclear matter at our calculated 
saturation density (binding energy) is about —18 MeV which is smaller than the empirical 
value (-16 MeV). 

For the polarized hot asymmetric nuclear matter, the nuclear symmetry energy is given 

by 

Fsym{p. r, 6n, 5p) = F(p, T, /? = 1, 5„, 6p) - F(p, r, /3 = 0, 5„, 6p). (24) 
In Fig. [3l we have plotted the nuclear symmetry energy {Fgym) of spin polarized hot asym- 
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metric nuclear matter as a function of the temperature at p = 0.31 and 0.16 fm~^ for two 
cases with Sn = Sp = 0.0 and Sn = Sp = 1.0. It is seen that at low temperatures, for 
a given polarization, the difference between the symmetry energies of different densities is 
substantially large. However, this difference decreases as the temperature increases. Fig. [3] 
shows that for each density, the symmetry energy increases by increasing the polarization 
and temperature. However, for each density, the difference between the symmetry ener- 
gies of different polarization decreases by increasing the temperature. From Fig. |3l for 
the case of unpolarized nuclear matter (5„ = 5p = 0), at T = MeV, we have found 
that the value of the symmetry energy corresponding to our calculated saturation density 
(p = 0.31 /m~^) is about 39 MeV. For this case of nuclear matter, the empirical value 
of the nuclear symmetry energy is 28 — 32 MeV corresponding to the empirical saturation 
density of about p = 0.16 fm~^ js^]- This shows that our calculated symmetry energy as 
well as our saturation density are greater than the empirical results. We have also found that 
the symmetry energy of the unpolarized nuclear matter at T = MeV for p = 0.16 fm"^ 
is about 26 MeV, which is less than its empirical value. The density dependence of the 
nuclear symmetry energy of polarized hot asymmetric nuclear matter has been also shown 
in Fig. m For each value of the polarization, we see that -Fsym is an increasing function 
of the density. However, the rate of increasing of the symmetry energy versus the density 
increases by increasing the polarization. 

The response of a system to the magnetic field is characterized by the magnetic suscep- 
tibility, X) which in the case of asymmetric nuclear matter is defined by a 2 x 2 matrix as 
follows, 



In the above equation, p is the total density, pi (pj) and pi (pj) are the magnetic moment and 
density of the particle i (j), respectively, and 6i and 6j are the spin asymmetry parameters of 
particles i and j, respectively. The onset of the spin instability of the system appears when 




(25) 



where the matrix elements 1/Xij ^I's given by. 




(26) 
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the sign of the determinant of this matrix becomes negative. We calculate the magnetic sus- 
ceptibility of the polarized asymmetric nuclear matter in terms of the ratio -^^j^^^-, where 
Xf is the magnetic susceptibility for an ideal Fermi gas containing noninteracting protons 
and neutrons. For the polarized hot asymmetric nuclear matter, the ratio dct(i/xF) ^^^^^^ 
temperature, for p = 0.16/m~^ at different values of the isospin asymmetry parameter 
has been presented in Fig. |5l This figure shows that for all values of the isospin asymmetry 
parameter and temperature, dtt(i/xF) positive. This behavior indicates that there is no 
magnetic instability for the hot asymmetric nuclear matter. From Fig. |5l we see that for all 
isospin asymmetry parameters, the ratio ^^^"^(x/xf) "^^ich is greater than unity, monotonically 
decreases by increasing the temperature. This shows that as the temperature increases, the 
strong correlation in nucleonic matter which arises by nucleon-nucleon interaction, becomes 
less important. Fig. [5] also shows that for all given temperatures, the ratio dtr(i/xF) 
creases by decreasing the asymmetry parameter This indicates that the proton fraction 
substantially affects the magnetic susceptibility of asymmetric nuclear matter. This effect 



has been already pointed out by Kutschera et al. [58[ and Bernardos et al. [20[. In Fig 



ini we have plotted the ratio dtt(i/x'^\ ^ function of the density at T = 20 MeV for a 
wide range of the isospin asymmetry parameters. We see that for all densities, this ratio 
is always positive, indicating that the spontaneous ferromagnetic phase transition does not 
occur for hot asymmetric nuclear matter at any densities. Fig. [6] also shows that this ratio 
is always greater than unity and always increases by increasing the density. Therefore, we 
can conclude that the strong nucleon-nucleon correlation becomes more important as the 
density of system increases. 

The equation of state of the polarized hot asymmetric nuclear matter for different values 
of the polarization and isospin asymmetry parameter at T = 10 and 20 MeV has been 
presented in Fig. [71 We see that for all temperatures and isospin asymmetry parameters, 
the equation of state becomes stiffer as the polarization increases. We also see that for 
a given polarization at each temperature, the pressure of spin polarized hot asymmetric 
nuclear matter increases by increasing the isospin asymmetry parameter. By comparing the 
two panels of Fig. [TJ we can see that increasing the temperature leads to the stiffer equation 
of state for spin polarized asymmetric nuclear matter. 
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IV. SUMMARY AND CONCLUSIONS 



We have used the lowest order constrained variational (LOCV) method to calculate the 
free energy of the polarized hot asymmetric nuclear matter at different temperatures, for 
various values of the isospin asymmetry parameter and spin polarization employing the AVis 
two-nucleon potential. Our results show that the free energy of this system increases by in- 
creasing both the isospin asymmetry parameter and the spin polarization while it decreases 
by increasing the temperature. We have also calculated the symmetry energy of this system 
to show that this term depends on the spin polarization, temperature and density. The mag- 
netic susceptibility of spin polarized hot asymmetric nuclear matter has been also calculated 
for a wide range of densities and isospin asymmetry parameters at different temperatures. 
Our calculations do not show any spontaneous phase transition to the ferromagnetic state. 
We have also computed the pressure of spin polarized hot asymmetric nuclear matter for 
different values of the polarization, isospin asymmetry parameter and temperature. Our 
results show that for all temperatures and densities, the equation of state of this system is 
an increasing function of spin polarization and the isospin asymmetry parameter. Finally, 
an agreement is seen between our results and those of other many-body calculations. 
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FIG. 1: The free energy per nucleon of spin polarized hot asymmetric nuclear matter as a function 
of the total number density (p) for unpolarized and fully polarized matter at T = 10 and 20 MeV 
for p = 0.0 (a), 0.3 (b), 0.6 (c) and 1.0 (d). 



14 




FIG. 2: The free energy per nucleon of the unpolarized hot symmetric nuclear matter (/3 = 0.0) 
versus the total number density {p) at T = 0, 10 and 20 MeV. 




FIG. 3: The symmetry energy of the spin polarized hot asymmetric nuclear matter versus temper- 
ature for two values of the density at different polarizations. 
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FIG. 4: The symmetry energy of spin polarized hot asymmetric nuclear matter versus density at 
T = 20 MeV for different polarizations. 
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FIG. 5: The magnetic susceptibility of hot asymmetric nuclear matter versus temperature at 
p = 0.16/m~^ for different values of the asymmetry parameter . 
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FIG. 6: The magnetic susceptibility of hot asymmetric nuclear matter versus the total number 
density (p) at T = 20 MeV for different values of the asymmetry parameter (P). 
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FIG. 7: The equation of state of spin polarized hot asymmetric nuclear matter for different values 

of the polarization and asymmetry parameter at T = 10 MeV (a) and T = 20 MeV (b). 
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